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Abstract 

We calculate threshold corrections to the running gauge and Yukawa couplings 
in the Exceptional Supersymmetric Standard Model (EgSSM) and analyse the 
more precise and reliable mass spectra in a constrained model (CEgSSM). Full 
expressions for the corrections are provided and the implementation into a 
spectrum generator is described. We find a dramatic reduction in the match- 
ing scale dependency of the masses of many states and observe a significant 
adjustment of the correlation of low-scale physical masses and high-scale pa- 
rameters. Still, in substantial regions of parameter space the mass of the 
lightest Higgs is compatible with the new boson discovered at the LHC and 
the model satisfies limits from collider searches for squark, gluinos and Z' 
bosons. We study the implications for gauge coupling unification from a new 
dependency of the spectrum on so-called survival Higgs fields which cannot 
be addressed without the inclusion of the threshold corrections. 
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1 Introduction 



Supersymmetric (SUSY) extensions of the Standard Model (SM) with TeV-scale 
SUSY breaking provide very attractive models for new physics which could be dis- 
covered at the Large Hadron Collider (LHC). Such models can solve the hierarchy 
problem of the SM due to the cancellation of quadratic divergences, enable embed- 
ding into Grand Unified Theories (GUTs), and thereby provide an explanation of 
the U(1)y rational charges, postulated ad hoc in the SM. These fundamental mo- 
tivations do not imply any minimality condition on the particle content or gauge 
structure. 

Here we consider the Exceptional Supersymmetric Standard Modef] (EgSSM) 
J2J [3], which is a non-minimal SUSY extension of the SM with an extra U(1)n 
gauge symmetry and new exotic matter at the TeV scale. A new Higgs singlet to- 
gether with an extra U(l) gauge symmetry solve the well-known /i-problem [1] of 
the Minimal Supersymmetric Standard Model (MSSM) without introducing the tad- 
pole or domain wall problems of the Next-to-MSSM^ [6j. These features also allow 
the lightest Higgs boson to be substantially heavier than in the MSSM or NMSSM 
already at tree-level, through new U(1)n -D-term contributions and an additional 
F-term contribution from the superpotential coupling between the singlet and dou- 
blet Higgs bosons, which can be substantially larger than the parallel NMSSM term 
due to different perturbative limits on the coupling [3]. 

The model is inspired and motivated by E§ GUTs, as the extra U(l) appears from 
the breakdown of Eq and the exotic new matter comes from complete Eq multiplets 
surviving to low energies, ensuring the cancellation of gauge anomalies. The extra 
£7(1) gauge group of the E 6 SSM is uniquely chosen such that right-handed neutrinos 
are neutral, allowing large Majorana masses and a high-scale see-saw mechanism. 

The Eq symmetry itself can arise from E$ x E' s Heterotic string theory after 
the breakdown of E's [7]. The E' s then plays the role of a hidden sector interacting 
with the visible sector only through gravitational interactions and leads to TeV-scale 
SUSY breaking, generating a set of soft SUSY-breaking parameters. 

In recent publications [H |9] a constrained version (CEgSSM) was introduced 
where all the soft SUSY-breaking masses are determined by a universal scalar mass, 
m , gaugino mass M1/2 and soft trilinear scalar coupling A. Owing to the unification 
of matter and Higgs fields in complete GUT multiplets the CE 6 SSM is particularly 



-'Tor a brief review see pQ. 
2 Reviews are given in [5]. 
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well motivated. In Refs. [8j [9] the low-energy mass spectra were explored for the 
first time with benchmark scenarios, representing the phenomenologically distinct 
possibilities in the model. This work has recently been updated to look at spectra 
which are consistent with the 125 GeV Higgs signal and recent searches for squarks 
and gluinos [10]. Other aspects of this model and similar variants have also been 
discussed in Refs. 11111121113111111511151113112111311211121]. 

The calculations in Refs. [8j |9j [10] neglect important threshold corrections and 
involve a significant dependence on unphysical threshold scales. This limits the 
accuracy of phenomenological results and it was stated that the precision for the 
masses was not better than 10%. 

However to facilitate model discrimination from data, and ultimately to recon- 
struct parameters in the event of a signal, precise predictions of the TeV scale SUSY 
masses predicted from GUT scale parameters are important. Indeed for the MSSM 
there are already a number of state of the art spectrum generators publicly avail- 
able [22j [23j EU [25j [26] which calculate the spectrum for scenarios, such as the con- 
strained MSSM (CMSSM). Typically these employ two- loop renormalization group 
equations, full one-loop matching conditions for gauge and Yukawa couplings and 
one-loop shifts to pole masses and comparisons of their predicted masses suggests 
an accuracy of about 1% [27]. 

The present paper is devoted to the calculation of threshold corrections to the 
running gauge and Yukawa couplings in the Exceptional Supersymmetric Standard 
Model (E 6 SSM), which enable a more precise evaluation of the mass spectrum from 
high-scale assumptions, and to an extensive study of the phenomenological conse- 
quences in the CEgSSM. As in Refs. [U [9j [10] we employ two-loop renormalization 
group equations for the running of the gauge and Yukawa couplings between the 
weak scale and the GUT scale. As a major improvement, we also compute and take 
into account the corresponding one-loop threshold corrections arising in the tran- 
sition from the full CEgSSM to the SM as a low-energy effective theory. We show 
that the individual threshold corrections indeed drastically reduce the dependence 
of our results on the unphysical threshold scale. The resulting, more accurate and 
reliable mass spectra are discussed. Particular attention is paid to the prediction for 
the Higgs mass and the gluino mass. We confirm the recent result of Ref. [10] that 
the lightest Higgs mass can easily be in agreement with recent LHC discovery of a 
new boson 2EJ EH]- Finally we use our improved precision to study gauge coupling 
unificatioro for the CE 6 SSM and the impact of the survival Higgs sector, which 



3 Previous studies (301 131] used trivial matching conditions, equivalent to assuming all new 
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appears only via threshold corrections. 

This paper is structured as follows. In Sec. [2Jwe briefly outline the E 6 SSM model, 
in Sec. [3] we describe the improved spectrum calculations, give analytical results for 
the threshold corrections. Sec. H] describes the application to an improved CEgSSM 
spectrum generator and explains the full numerical procedure used. In Sec. |5] we 
illustrate the effect of the thresholds and the improvement in accuracy of our results, 
and we give an extensive discussion of the resulting model predictions. 

2 The E 6 SSM 

The E 6 SSM is a supersymmetric gauge theory, inspired by Eq GUTs in which the 
Eq is broken at the GUT scale Mx via the Hosotani mechanism [32] down to the 
low-energy gauge group of the E 6 SSM, 

SU(3) C x SU(2) L x U{l) Y x U(1) N . (1) 

Here the £/(1)jv is a special case of a U(l)' symmetry arising from Eq breaking, 

Ef(l)' = tf(l) x cos0 + i7(l)^sin0, (2) 

with tan# = vT5, where the gauge groups U(l) x and U(l)^ are defined by the 
breaking of Eq -> 50(10) x U(l)^ and 50(10) ->■ SU(5) x U(l) x [331 132] • 

The choice of £7(1) at makes the right-handed neutrino a pure gauge singlet, thus 
allowing a gauge invariant Majorana mass term enabling a high-scale seesaw mecha- 
nism for the generation of neutrino masses. All other matter from complete Eq mat- 
ter multiplets survives to low energies, thereby ensuring the cancellation of gauge 
anomalies. 

The three families of E 6 SSM matter particles fill complete (27); representations 
of the E 6 group, ensuring that full low energy gauge group, including the new U (1) ^ 
gauge symmetry, is anomaly free. In addition to these, the model has two Higgs- 
like doublets H' and H', the so-called survival Higgs doublets, both originating 
from extra (27)' and (27)' representations to ensure gauge coupling unification at 
a high scale Mx- The decomposition of the fundamental (27) representation under 
SU(5) x U(1)n and (JTJ is listed in Tab. [TJ All Standard Model fermions as 
well as their superpartners fit into the multiplets (10, 1); and (5, 2)j. The (5, —3) 
and (5, —2) representations contain Higgs-like doublets Hu, H 2 i and exotic colored 

EgSSM states and all new MSSM states have common masses, Te 6 ssm and Tmssm respectively. 
Such scenarios are clearly not realised in the constrained version of the model. 
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Tab. 1: EgSSM SUSY multiplets and their gauge quantum numbers (gen- 
eration index i = 1,2,3), where G S M = SU(3) C x SU(2) L x U(1) Y - For the 
abelian groups U(l)y and U(1)n the charges Y/2 and N/2 are listed. 



matter fields -D,, D,. The remaining SU(5) singlets (1,0), and (1,5), equate to 
right-handed neutrinos and fields S,, respectively. A complete particle listing can 
be found in Tab. [2J 

In an Eq GUT the gauge bosons and their superpartners fit into the adjoint (78) 
representation of the E e , which is decomposed under the low-energy gauge group 
SU(3) C x SU(2) L x U(1) Y x U(1) N into 

(78) (8, 1, 0, 0) + (1, 3, 0, 0) + (1, 1, 0, 0) + (1, 1, 0, 0) + • • • (3) 

The gluons present in a low-energy model are associated to (8, 1, 0, 0). The (1, 3, 0, 0) 
multiplet contains the weak gauge bosons and the two U(l) gauge fields belong to 
the two (1, 1,0,0) representations. When the E$ is broken at the GUT scale Mx, 
the other gauge bosons are expected to have masses of the order of Mx- However, 
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Tab. 2: Component fields of the EgSSM superfields (generation index % = 


1, 2, 3). The charge conjugation of a spinor is defined as 


:= C ip T , where 



C = i 7 2 7°. 

the E 6 SSM is a low-energy model by construction and does not include these heavy 
GUT gauge bosons. 

As with the MSSM embedded into GUT models, the new bosons absent from 
the low-energy theory could give rise to significant threshold corrections to gauge 
coupling unification at Mx- This should be borne in mind when studying the success 
or failure of gauge coupling unification within a low-energy SUSY model. However 
the purpose of this paper is to consider the low-energy threshold effects from sparticle 
masses and to test the significance of these on the physical predictions made by 
postulates about the high scale parameters. 

For (27)j representations of E$ the most general renormalizable superpotential 
with full Eq invariance is given by the trace of the (27) i x (27) j x (27)k- Invariance 
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under the low-energy gauge group allows further terms. But as in the case of the 
MSSM, the most general gauge invariant superpotential is not phenomenologically 
viable [2J [8] , as it contains baryon number violation and unacceptably large contri- 
butions to non-diagonal flavor transitions. To conserve baryon number and avoid 
flavor changing neutral currents, one first imposes a Z^ symmetry, under which 
an « superfields transfer as odd, except H a , tf 23 and sM 

The remaining 

S77(3) c x SU(2)l x U(l)y x Z7 (1)jv an d invariant superpotential reads 

^e 6 ssm = -y^HrsL^ef - y*(H 13 Qi)$ - y^Q^uf 

+ \M l3 N?Nf + h%{H lz H>)e° + h%{H^H')Nf + fi'(H'H') (4) 

+ \S 3 (HuH 2 i) + KiS 3 DiDi + f a i3S a {H 13 H 2 p) + fapSjJtxpH^), 

where the SU(2) superfield spinor product is defined as (AB) := A 2 B 1 — A 1 B 2 . 
A problematic consequence of this Z^ symmetry would be that the exotic quarks 
would only have gauge interactions and interactions with the singlet field, leading to 
stable charged matter in violation of experimental constraints [35j [361 EZl • Therefore 
the Zif symmetry can only be approximate. 

The Z^-violating terms should not lead to rapid proton decay. Hence, another 
discrete symmetry, analogous to i?-parity, must be required. This can be done in 
two ways: either a Z% symmetry is implemented, under which all superfields except 
the leptons are even (Model I) or one imposes a Z^ symmetry, under which the 
exotic quarks and leptons are odd whereas the others remain even (Model II). The 
superpotential (J4]) is then enlarged by one of the following Z|^-violating but Z^ 3 - 
conserving terms 

WWi i = 9? jh bi{QjQk) + 9 q ijk bidfuk, (5) 
^Modoi ii = gfjkNfDjuC + gf^efbju^ + g? k (QiLj)b k . (6) 

In Model I the scalar exotic quarks can decay into two quarks (they are diquarks) 
and in Model II they are leptoquarks since they decay into a lepton and a quark. 

For correct electroweak symmetry breaking only the scalar components of H13, 
H23 and S3 get a non-zero VEV. To ensure this, we impose that a certain hierarchy 
between the Yukawa couplings must exist 

«t ~ A3 > Ai )2 > /a/3, /a/3, hfj, h^. (7) 

4 Although a high-scale family structure is not a part of our model by construction, such a 
symmetry could be the result of a A27 family symmetry at the GUT scale [15]. 
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This hierarchical structure allows a simplification of the superpotential (J4]). Inte- 
grating out the right-handed neutrinos, which are assumed to be very heavy, and 
keeping only the dominant terms, one arrives at 

^e 6 ssm ~ -yr(HdL 3 )e% - y b {H d Q 3 )d% - y t {Q 3 H u )u 3 

(8) 

+ XiSiHuHti) + KtSDiDi + fJ.'(H'H'), 

where the scalar components of H u := H 23 and H d := Hi 3 and S := S 3 develop 
VEVs, (H®) = v u /y2, (H%) = Vd/y/2, giving mass to ordinary matter while (S) = 
s/y/2 gives exotic quark masses, KiS — > Kis/^/2 =: fj,^; masses for the fermion 
components of the first two generations of "inert" Higgs-like doublets (the ones which 
don't get VEVs) fig := A q s/a/2 and an effective //-term for the Higgs doublets 
/i c fj := \ 3 s/y/2. Eq. (jSJ) is the superpotential which will be used in the following 
analysis to determine the particle spectrum from high scale assumptions, inspired 
by minimal Supergravity and Eq GUTs. 

The Higgs potential, electroweak symmetry breaking conditions and mass eigen- 
states of all the particles in the model have been presented in Ref . [8] . 

3 Threshold corrections in the EgSSM 

The E 6 SSM is a low-energy model which is motivated by a particular high-scale 
structure. The high and low scales are connected by renormalization group equations 
which have already been given at the two- loop level in Ref. [8]. In the following we 
present the results for the threshold corrections which are required for a consistent 
inclusion of subleading effects. Before presenting our computation and the results 
we give a brief summary of relevant general features of threshold corrections. 

3.1 Threshold corrections for matching full and effective 
theories 

To connect running fundamental EgSSM parameters with SM quantities we consider 
the SM as an effective theory of the EgSSM and match the two theories at a threshold 
scale T E6 s SM , which should be of the order of the heavy E 6 SSM particles [5H] . To 
see the need for and the properties of threshold corrections most clearly, we slightly 
generalize and suppose we have a full and an effective gauge theory and we want to 
calculate the DR gauge coupling gMi(Q&x) at a fixed scale Q^ x in the full theory. As 
input we know g e s(Qiow) in the effective theory at some low scale Qi ow < Qfi x - We 
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start with our effective theory at Q\ ow and run the gauge coupling g eS from Q low to 
the matching scale Ti, using the beta function /3 e g of the effective theory, see Fig. [TJ 
At Ti we match the full and the effective theory and calculate gf u u(Ti) without the 
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Fig. 1: Visualization of gauge coupling threshold corrections 



use of threshold corrections, 



0fuii(7i) := £? c fr(Ti 



(9) 



Now we use the beta functions of the full theory (3 {u u to evolve <7f u n(Ti) to the desired 
scale Qflx and obtain gMi(Qnx)- 

The problem with this approach is that the resulting coupling <7fuii(Qfix) depends 
on the choice of the unphysical matching scale Ti, i.e. gf u \\(Q&x] Ti). More precisely, 
if one would do the same calculation using another matching scale T 2 the difference 
would be 
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because the solution of the RGE is given by 
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which means that the gauge coupling in the full theory at the scale depends 
on the unphysical matching scale t and the slope of g^(Qftx] t) with respect to t 
is proportional to the difference of the beta functions of the full and the effective 
theory. 

The way out is to use threshold corrections Ag when matching the full and the 
effective theory, i.e. replace Eq. by 

SfuiiC^) := ^(TO + A^Ti), (13) 

where Ag(Tj) can be obtained by matching Green functions of the full and the 
effective theory, see below. As a result the right-hand side of Eq. f[T2~j) will vanish, i.e. 
the gauge coupling <7f u ii(Qfix) w iU n °t depend on the matching scale. This property 



will be used as a test in Sec. 15.11 

To calculate the threshold corrections we consider a full theory consisting of 
parameters pi, . . . , p p , light fields l\, . . . ,l q and heavy fields Hi, . . . , H r 

£ full = £ full (pi) ; pp . ^ ^ ^ jj u ; Rr) (14) 

The effective theory might contain effective parameters p±, . . . , pk and effective light 
fields li, . . . , lg, where the heavy fields were integrated out 

C cS = C cS ( Pl ,...,p k ;h,...,i g ). (15) 

The requirement that both the full and the effective theory describe the same physics 
in the limit p — > can be achieved by equating all connected Green functions with 
light external fields li in the full and in the effective theory at zero external momenta. 
This condition leads to the equality of all one-particle irreducible correlation func- 
tions T that are one-particle irreducible with respect to the light fields U (1LPI) 
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and inserting Eq. ([IT]) - (119"]) into (fT6|) yields the matching condition 

pfulltree pJull.lL = f j + * y- ^ ] peff .tree + pefflL ^ + (2Q) 

This equation is evaluated at zero external momenta. The analogous equation holds 
for the derivatives of the 1LPI correlation functions. Imposing all matching condi- 
tions yields the definitions for the K[. in terms of 1LPI correlation functions and at 
the same time the desired relations between the parameters of the effective and full 
theory (threshold corrections) 

Pl = Pi + Ap 4 (1 - loop) ( Pl , . . . , Pp , K h , . . . , K lq ), (i = 1, . . . , k). (21) 



3.2 Gauge coupling threshold corrections in the E 6 SSM 

Strong gauge coupling To calculate the threshold correction for the gauge cou- 
pling g 3 of the unbroken SU(3) C one needs to apply the matching procedure of 
Sec. 13. II to the following 1LPI correlation functions 

d P Itf M fe "P)U = 9 ? r M>' -*)U > ( 22 ) 
9* r aSi^ - k ) , 2 _ n = ^ T ^Gi(k, -k)\ k2=Q , (23) 
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pSM 

p=fc=0 

where are the light colored fields that remain in the SM, i.e. the SM quarks. 
After imposing the relative field renormalizations (|T9|) and decomposing the 1LPI 
functions into a tree-level and a loop part one gets the threshold correction for g 3 

g msu = ^ms,e 6 ssm + _|_ ^ * FfC (r f ) log + \F.C{r.) log ^ 

(25) 

The sums run over all heavy fermions / and scalars s that are integrated out. 
Eq. ([25]) is also in agreement with the general result in [39]. The constants C(r) are 
invariants of the representations r of SU(3) C and are given by 

C(N) = - (fundamental representation N), (26) 

C(G) = 3 (adjoint representation G), (27) 

and Ff, F s account for the different number of field degrees of freedom 

1 if / is a Dirac fermion, 
1/2 if / is a Majorana fermion, 
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(28) 
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if s is a complex scalar, 
1/2 if s is a real scalar. 



(29) 



When matching the E 6 SSM to the SM, only the gluino, the squarks and the exotics 
contribute in Eq. fl25f) and we have 
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where we have added a finite counterterm which converts #3 from the MS to the DR 
scheme 



Electroweak gauge couplings The calculation of threshold corrections for the 
couplings qy and g 2 is more involved, because the gauge group SU(2) L x U(l)y is 
spontaneously broken to U(l) em . From the relations 
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one can see that the threshold corrections for gy and g 2 are related to those of the 
W ± and Z boson masses as well as to the gauge coupling e of the remaining U(l) em 
gauge symmetry. Therefore the following matching conditions are imposed to obtain 
threshold corrections for e, mz and mw 
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The additional matching condition for T Zfl A„ is necessary, because the gauge fields 
By, mix to Ap, and W^. Introducing relative field renormalizations for W^, 
Zp, Ap, ipiL and ipm and inserting them into the matching conditions fl34"l) - fj4"0]) leads 
to the threshold corrections for the W and Z boson masses and the electromagnetic 
coupling^ 



m 1 



SM\2 
V I 



[m 



E 6 SSM\2 

V J 

E 6 SSM\2 

V ) 



+ r 



E 6 SSM,lL,heavy 
VuVu,T 



fc 2 =0 



- myK vv , 



Am 



Vi 



3 E 6 SSM 



-K 



s w 



AA 



K 



ZA 



2 2c w 

The relative field renormalization constants in Eq. ( 141 j) and 



Ve{W,Z}, (41) 
(42) 

E 6 SSM + Ae _ (43) 



are given by 



Kaa 
K ww 







^E 6 SSM,lL,heavy 



Q k 2 A»A„,T 



k 2 =0 







^E 6 SSM,lL,heavy 



dk 2 W?W^,T 



K ZA 

Kzz 



m 



pE 6 SSM,lL,heavy 
2 1 A„Z„,T 



9 



fc 2 =0 
E 6 SSM,lL,heavy 



ZuZv,T 



fc 2 =0 



(44) 
(45) 



fc 2 =0 ' <9/c 2 
and by r^eSSM.iL.heavy we denote the one-loop parts of the gauge boson 2-point 
functions that contain all particles that we integrate out, i.e. all non-SM particles. 
Using Eq. (|4T]) -( l43|) one can write the threshold corrections for gy and g% in the 
form 



9y M 



9^ 



J3M 



//l7 



SM 
W 



in 



EgSSM 

9y 



e 



1 Am| 1 Am 2 ^ 



2 m| 



2 m 2 ^ 



,SM 



EeSSM 
i/2 



Ae _ c 2 w Am| 



+ 



2m 2 



(46) 
(47) 



Inserting the explicit form of Ae, Am^ and Am| yields 

3 



9 Y 



DR,E 6 SSM 



MS,SM 



a? 




8=1 



5 The Ward identity which reflects the U(l) em gauge invariance in the Standard Model was used 
to simplify the result. 
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/ J- i 



^EEE 



m 



log 



i=l p=l j=l 
9 3 2 2 /y 

;EEEH 



/' 



3=1 



#23 \ i m # 



log 



1=1 p=l j = l 




-5E 



3=1 



3 ^ \ 2 

3=1 V 



^ d i=i 



2 /v x 2 _ r, 2 /y =J \ " m = 



E 




log 



log 



if 



(48) 



DR,E 6 SSM MS,SM . #2 
92 = 92 + ~71~^ 

1 2 (4vr) 



5-^E>^-;E> 



1=1 



4. 

log — — 

3 yU 6 



u 6 z — ' 

^ i=i 

1 m# 



2 2 



EE'°< 

i=i p=i 



1 



3 2 



aEE'o^ 

8=1 p=l ^ 



1 , m H> 



1 m H / 

- log log 

pL 3 /X 

1 mj 1 

- log — - - log ^ 

O /i 3 /! 



(49) 



where we have added a finite counterterm which converts g% from the MS to the DR 
scheme [40] and neglected mixing. 

3.3 Yukawa couplings in the E 6 SSM 

Since the SM fermion masses are measured but not the Yukawa couplings, we don't 
define the EgSSM Yukawa couplings in terms of threshold corrections to the running 
SM Yukawa couplings. Instead we directly match them at the one-loop level to the 
measured SM fermion masses and gauge couplings via 



Vt 



DR,E 6 SSM 



92 



DR,E 6 SSM on-shell,SM 



m 



5m 



m 



on-shell,SM 
W 



sin (3 



on-shelLSM 



+ 



5m t 



m 



on-shclLSM 



(50) 



DR,E 6 SSM 



Vb 



DR,E 6 SSM 



92 



DR,E 6 SSM DR,(5) 



5m 



m 



on-shcll,SM 
W 



cos (3 



on-shell,SM 



+ 



Srrih — 5m 



(5) 



5m shift 



m 



DR,(5) 



(m z ) 



DR,E 6 SSM on-shell,SM 

92 'J^t 

[K on-shcll,SM n 

\J 2m w cosp 



5m 



on-shell,SM 



+ 



5m 7 



ni 



on-shell,SM 



(51) 
(52) 
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Here m ° n " shell ' SM ; m on-sheii,SM anc j m ^~ shell ' SM are Standard Model on-shell masses 



To avoid large logarithms from the bottom mass we use the DR value of mj in the 
5-flavor QCD at mz, rn^ K (mz) = 2.83 GeV [12] which we then shift to the scale 
ji where the Yukawa couplings are evaluated at. The counterterms in Eqs. fl50l)- (l52l) 
are defined as 



(f = t,b,T), (53) 

finite 



A DR,(5)/ \ / 2 

Sm f = K «-.''.QCD (5) _ fn p m = J_^_M gl ( 5 + 3 log ^ ) , (54) 



,2 



Sm s r = ^,(5) lQg ±_ = _2m^ 5 \m z )fL l og JL (55) 

b mz 4:ir z mz 

r r on-sheIl,E 6 SSM r DR,E 6 SSM rr / 2 \ /r R \ 

om w = om w —dm w ' = MeYi ww ^ T {m w ) . (56) 

finite 

The self-energies and Ily^jy.T are listed in Appendix O It was checked that the 
divergences of (I50l) - (l52|) . see Appendix [Dj are in agreement with the prediction from 
the one-loop RGEs. 



4 Obtaining accurate spectra in the CEgSSM 
4.1 The CE 6 SSM and its parameters 

The threshold corrections presented in the previous section improve the precision of 
the high-scale-low-scale connection in the general EgSSM. In the present paper we 
apply them to the constrained version of the E 6 SSM (CE 6 SSM). The CE 6 SSM is 
defined by a universal scalar mass m , gaugino mass M\n and trilinear mass A at 
the gauge coupling unification scale Mx E] . Owing to the unification of matter 
and Higgs fields in complete GUT multiplets these constraints are particularly well 
motivated. The soft scalar masses from the survival Higgs (which appear in the 
incomplete (27)' and (27)' representations) sector, m H >, mg, are not assumed to be 
unified with mo and the soft bilinear Bfi' is also unconstrained. 

The CEgSSM is an example of a highly predictive, non-minimal SUSY model. 
It contains many new states at the TeV scale but has only few free parameters. 
In theory, the CE 6 SSM is fully determined by specifying the gauge couplings, the 
superpotential parameters yt,b,r, K and K\ and the universal soft parameters mo, 
Mi/2, A at the GUT scale as well as the survival Higgs parameters. All high-scale 
parameters are shown in the top row of boxes in Fig. [21 

In practice, the model predictions of course need to agree with known SM con- 
straints, in particular with the four known masses mz, th,6,t and the known low- 
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energy gauge couplings. Furthermore, it is useful to take the low-energy values of 
tan/3 = v u /vd and s = y/2(S) as free parameters of the model. In this way, the 
six GUT-scale parameters mo, M1/2, A and yt,b,r are traded for the four known SM 
mass constraints and the two low-scale input parameters tan f3 and s. Compared to 
the familiar CMSSM, the CE 6 SSM has no fi and B[i parameters, which could be 
adjusted to fulfill SM constraints. Hence, the CMSSM input parameters m , M 1 / 2 
and A are calculable in the CEgSSM. 

Fig. [2] shows the resulting structure of input and output of the CE 6 SSM and the 
connection of the high and low scales. High-scale input parameters are the super- 
potential parameters Xi(M x ) and Ki(M x ); low-scale input parameters are s(Qfi x ), 
tan /3(Qfi x ) and the survival Higgs parameters (defined at the matching scale Te 6 ssm), 



all defined in the DR scheme. The low-scale SM constraints fix the remaining pa- 
rameters. In previous studies Aj, Kj, s and tan/3 were sufficient to fix the spectrum. 
Here, with the inclusion of threshold corrections, we will be able to study the per- 
turbation caused by the survival Higgs sector as well. 

As the Figure shows, 6 + 3 + 4 = 13 high-scale parameters have been traded 
against these 13 low-scale input parameters, so this high- and low-scale input de- 
termines the structure of the model completely and at all scales. The fixed scale 
Q fix , where s and tan/3 are defined at, is set in Sec. \5\to a value of the order of the 
matching scale Te 6 ssm- 



high-scale input: 



low-scale input: 

„SM ™SM 

9i ! ! 

m z , t/3, s, 
fi', m H >, m s ,, Bfj! 



unification scale M 



x 



E fi SSM RGEs 



matching scale T\ 



E,iSSM 



5 f« SSM ,yf SSM = function of SM 
parameters and E 6 SSM spectrum 



high-scale output: 

E 6 SSM E 6 SSM 

Hi j y/ j 
m , My s , A, 

fi', m H ', rag,, Bfi' 



low-scale output: 
A;, Ki, £ soft , spectrum 



Fig. 2: The structure of high- and low-scale input and output of the CEgSSM. 
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4.2 Numerical procedure and the improved spectrum gen- 
erator 

We extended the particle spectrum generator previously written for [8], which is 
partly based on SOFTSUSY 2.0.5 [21]. As input the program gets a CEgSSM param- 
eter point specified by the high- and low-scale input 

Ai,2,3(M x ),K li2i3 (M x ),s(Q fix ),tan/3(g fix ), 

(57) 

H {TE 6 ssM),mH'(T E6SSM ),mfi,(T E6SSM ), B/i (T Ef , SSM ), 

and the program knows about the SM constraints discussed above. As illustrated in 
Fig. [21 RGEs and threshold corrections connect high and low scales, and the aim is 
to find output parameters consistent with the given input and then to compute the 
physical particle masses from the obtained low-energy SUSY breaking parameters. 
However, the threshold corrections can only be computed once the full mass spec- 
trum is known. Hence, these steps need to be iterated until convergence is reached, 
and in the first iteration the threshold corrections must be ignored. 

In the actual computation, we divide the basic strategy of each iteration into 
five steps, see Fig. El The details are as follows. 





determine 

„E 6 SSM „E 6 SSM . T 

fM ,Vf at i Ee SSM 








SUSY iteration 








determine soft pa- 
rameter dependency 








EWSB iteration 








mass spectrum 












threshold corrections 





Fig. 3: Program flow chart. The dashed box marks the structure of the old 
spectrum generator, which had no threshold corrections. 

1. Determine the gauge and Yukawa couplings in the EgSSM at the threshold 
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scale T E6 gg M from the known SM gauge couplings and masses, using SM RGEs 
and the threshold corrections. 



(a) Evolve the SM MS gauge and Yukawa couplings from their known values 
at m z [H] to the intermediate matching scale T E6S g M using SM RGEs. 



(b) Convert the SM MS couplings to E 6 SSM DR couplings using the indi- 
vidual particle threshold corrections discussed in Sec. |3j This step is 
only possible once the physical particle masses are known, i.e. after the 
first iteration. Hence, in the first iteration, the threshold corrections are 
replaced by the trivial definitions 

„sm SM 

yf SSM --=—Af = b > t)> y? 6SSM ■= ^ ( 58 ) 

J cos p sin p 

? WM :=fl p M(i = 1)2)3)i (59) 

Note that, in contrast to Ref. [8], we do not use an intermediate matching 
to the MSSM because the MSSM and the E 6 SSM spectra are typically 
mixed. 

2. Use RGEs to determine the high-scale gauge and Yukawa couplings and the 
low-scale values of the Aj, «j. This step is completely independent of all soft 
SUSY breaking parameters because of the structure of the RGEs. 

(a) Estimate low-energy values for the Yukawa couplings Aj (T Ee . SSM ) , «j (Te 6 ssm) 
The necessity of this estimate is an additional complication not present 
e.g. in the CMSSM. The difference is due to the fact that in the CMSSM 
all high-scale input parameters are soft breaking parameters, while here 
they are superpotential parameters. 

(b) Evolve all E 6 SSM gauge and Yukawa couplings from T EfiSSM up to the 
unification scale Mx, defined as the scale where g± — gi, using two-loop 
E 6 SSM RGEs. 

(c) Set Yukawa couplings \(Mx), K>i(Mx) to program input values. 

(d) Set := g /\/40, where go is defined by ji(Mj). 

(e) Perform an iteration between Mx and Te 6 ssm to obtain values for the 
gauge and Yukawa couplings and Aj, Kj until consistency is reached with 
low-energy boundary conditions yf(T E( . SSM ), gi(T E( . SSM ) and high energy 
boundary conditions Aj(Mx), Ki(Mx)- 
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3. With the gauge, Yukawa and X i: couplings now known at all scales we 
find solutions for the soft breaking parameters. They are determined by the 
high-scale constraints of universal mo, M1/2, A, and by low-scale electroweak 
symmetry breaking (EWSB) conditions, i.e. consistency with the measured 
value of rriz and the input values for tan/3 and s. First, the dependency of 
the low-energy soft mass parameters on the GUT scale values wip, Mu 2 , A is 
expressed in terms of the semi-analytical formulas 

m^t) = Oi(t)ml + k{t)Ml /2 + Cl (t)A M 1/2 + (U(t)J%, (60) 
Ai{t) = ei {t)A + fi{t)M 1/2 , (61) 
Mi(t) = Pi (t)A + qi{t)M lf2 , (62) 

where the coefficients a>i(t) , . . . , qi(t) are calculated numerically at the scale 
i = logT E6SSM /M x . 

4. Next, the EWSB conditions are used to fix the values of the soft parameters. 

(a) The obtained expressions for the low-energy soft parameters mf (t), Mj(t), 
Ai(t) are then combined with tree level EWSB conditions 

dV_ _ dV_ _ dV 

dvi dv 2 ds 

and the known values of mz, tan/3 and s to form three quadratic equa- 
tions in the soft masses mo, Mi/ 2 and A, which can be reduced to a single 
quartic equation. 

(b) This quartic equation is then solved numerically and the values used to 
determine the mass spectra. Note that in principle there will be four 
solutions, though some or all may be complex. Therefore our routine 
deals with between zero and four sets of real solutions. 

(c) Tadpole corrections can now be calculated and added to the EWSB condi- 
tions, and a solution of consistent EWSB for the leading one-loop effective 
potential is found iteratively. 

5. The full CEgSSM particle mass spectrum is now determined for each set of 
{toq, M1/2, Aq} found to be consistent with both high scale and low scale 
boundary conditions. 

The final solution is found by iterating over all five steps until convergent solu- 
tions are obtained. 
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It is important to note that in general EWSB is not guaranteed in the CE 6 SSM, 
i.e. solutions for {wq, Mi/ 2 , A } from Eq. are not always found. But for suffi- 
ciently large values of Ki the soft parameter m| always gets negative at low energies, 
which triggers EWSB [8]. 

5 Results 

In this section we present the results obtained with our improved CE 6 SSM spectrum 
generator. We start by quantifying how the threshold corrections stabilize the results 
for parameters like gauge couplings and the mass spectrum. We will then extensively 
discuss the resulting, more accurate mass spectrum. On the one hand we will obtain 
more precise information on previously defined benchmark points, implying that 
some are now excluded by LHC data. On the other hand we scan over the parameter 
space and find large regions which are compatible with LHC limits on SUSY particles 
and are consistent with the latest discovery of the new boson which we associate 
with the lightest Higgs boson in the model. Finally we investigate the impact of 
the survival Higgs sector, which is possible due to the threshold corrections. It has 
implications on gauge coupling unification and on the predictions for the low-energy 
mass spectrum. In subsequent studies we sometimes make use of a test point, PP1, 
to illustrate effects, which is defined by, 

tan/3(Q fix ) = 35, A 1)2 , 3 (M X ) = k w {M x ) = 0.2, s(Q fix ) = 10 TeV, 

( 64 ) 

/A?e 6 ssm) = m H ,(T E6 ssM) = m s ,(T E6 ssM) = 10 TeV, -B//(Te 6 ssm) = 0. 
Furthermore we set = 3 TeV for all the following analyses. 

5.1 Matching scale dependency 

In Sec. 13.11 we presented Fig. [TJ to illustrate how the threshold corrections should 
adjust the RG flow of the gauge couplings so that the matching scale dependency 
is removed. Now in Fig. H] we present an explicit demonstration of this effect in the 
E 6 SSM with the threshold corrections we have calculated. The gauge couplings are 
plotted as a function of the renormalization scale for two different matching scales, 
T\ = 500 GeV and T 2 = 10 TeV. If no threshold corrections are used (upper plot) 
different matching scales yield different predictions of Mx, which leads to a different 
phenomenology at the TeV scale. The inclusion of threshold corrections reduces 
this unphysical behavior. As shown in the lower plot, it leads to an approximately 
matching-scale independent prediction of Mx- 
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This improved behavior should manifest itself in a reduced scale dependence of 
model predictions at the TeV scale. As a test of this, Fig. \5\ shows the E 6 SSM gauge 

without threshold corrections 

5.0 
1.5 
1.0 
3.5 
3.0 

-5 2.5 

i— I 

2.0 
1.5 
1.0 
0.5 
0.0 

10° 10 2 10 4 10 6 10 8 10 10 10 12 10 14 10 16 10 18 
renormalization scale / GeV 

inclusive threshold corrections 

5.0 
1.5 
1.0 
3.5 
3.0 

~5 2.5 

i — I 

2.0 
1.5 
1.0 
0.5 
0.0 

10° 10 2 10 4 10 6 10 8 10 10 10 12 10 14 10 16 10 18 
renormalization scale /i / GeV 

Fig. 4: Running gauge couplings for parameter point PP1 for two different 
matching scales T\ = 500 GeV and T2 = 10 TeV with and without threshold 
corrections. 

couplings gi(Q 'ax) at the fixed scale Qfi x = 3 TeV, as a function of the matching scale 
^e 6 ssm for parameter point PP1. In the case of the trivial matching relations ( 159|) one 
finds a clear (unphysical) dependence of QiiQ^) on T E6 ssm- As shown in Eq. ( 1T2"|) the 
slopes of (4vr) 2 /[2^(g fix; T E6SS m)] are given by the difference Aft := /3f eSSM - (3f M 
of the one-loop gauge coupling beta functions, up to two-loop and iteration effects. 
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Fig. 5: Dependency of the gauge couplings g. L at Q nx = 3 TeV on the matching 
scale Te 6 ssm for parameter point PP1. The circles show the behavior without 
threshold corrections and the squares with corrections for gi and yf. 



In Tab. [31 columns 4 and 5, the numerical values of A/3j as well as the slopes 
without threshold corrections for PP1 are listed. The values roughly coincide, thus 
confirming Eq. (TT2"j) . When the threshold corrections (T4"gj) . (I49p and (13 2 p are used, the 
matching scale dependency reduces to about 4 % or less, resulting in the flattening 
of the curves shown in Fig. |5j The remaining matching scale dependency as well 
as the difference between columns 4 and 5 are due to still missing contributions of 
higher orders. 
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Coupling 


^E 6 SSM 


/3f M 


Aft 


slopes of (4:7r 
w/o thresh. 


)7[2o } 2 (Q fix ;T Er gs M )l 
w/ thresh. 


9i 


9.6 


4.1 


5.5 


5.67 


0.18 


92 


4 


-3.17 


7.17 


7.42 


0.32 


93 





-7 


7 


7.33 


0.29 




oEeSSM frp \ 

Pf Ue 6 ssmJ 


/3/ M (^e 6 ssm) 


"A/3; 


slopes of (Att' 
w/o thresh. 


2 losruf (Qfi x ; TrrSSm) 
w/ thresh. 


Vt 


-2.33 


-6.14 


-3.77 


-3.81 


0.77 


Vb 


-4.66 


-8.18 


-4.92 


-4.95 


-0.20 


Vt 


-0.34 


0.77 


-0.29 


-0.21 


0.53 



Tab. 3: Effect of the threshold corrections on the dependency of the gauge 
and Yukawa couplings on Te 6 ssm for parameter point PP1. The slopes in the 
last two columns are obtained by linear fits to the data in Fig. [5] and El 



Fig. shows the same analysis for the Yukawa couplings y t , y^ and y T . Analogous 
to the gauge couplings one expects the following slopes of the Yukawa couplings as 
a function of the matching scale 

dt ~~ ~ (4tt)2 ^ - Pt W + /WW cos P\ = ~7^ A ^' 

(65) 

dlogZ/fe,r(Qfix;t) 1 foEeSSM/.N oSM/^N a ( + \ -2 ol _ 1 a q 

dt ~ ~ (4tt)2 L^V - P&,t I* J - /WW sm 0| = -— -A£ 6;T . 

(66) 



Note that these relations also contain the (5 function for tan /3, defined in Eq. f lD.14j) . 



In Tab. [3j columns 4 and 5, the numerical values of A/3j as well as the slopes without 
threshold corrections for PP1 are listed. One finds an agreement with the predicted 
slopes within a few percent. When the threshold corrections from Sec. 13.21 and 



13.31 are taken into account one finds an overall reduction of the matching scale 
dependency effects. The remaining scale dependence is due to neglected higher- 
order effects, including in particular QED-logarithms of the small fermion masses 
in higher powers than taken into account in the one-loop computation of 8m f in 
Eqs. ((5DD-(J52D. 
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Fig. 6: Dependency of the Yukawa couplings y/ at Qf\ x = 3 TeV on the 
matching scale Te 6 ssm for parameter point PP1. The circles show the behavior 
without threshold corrections and the squares with corrections for gi and yf. 



5.2 Particle masses 

In softly broken SUSY models, like the CE 6 SSM, where the pattern of soft breaking 
masses is set substantially above the EW scale (e.g. Mx), the low-energy DR running 
masses are very sensitive to the renormalization group flow of the gauge and Yukawa 
couplings. This implies that the physical mass spectrum strongly depends on the 
matching scale used for the gauge and Yukawa couplings unless stabilized by the 
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inclusion of threshold corrections. Here we show how the inclusion of the CE 6 SSM 
threshold corrections improves the prediction of the low-energy DR soft masses and 
yields a mass spectrum where the dependence on the unphysical matching scale is 
considerably reduced. 

In Figs. [7] and [8] the dependency of the particle masses on the matching scale for 
parameter point PP1 is shown. The matching scale is varied in the range [|T , 2T ], 
where To = 1.9 TeV is the geometric average of all particle masses shown in Fig. |HJ 

Fig. [7| focuses on the gluino. The theory uncertainty implied by the Te 6 ssm 
dependence is 65 % without threshold corrections (the percentage value is defined 
as the full variation of the particle mass divided by the mass at T Q ). This huge 
uncertainty is entirely due to the missing threshold corrections, and including these 
reduces the uncertainty to only 0.5 %. 

Fig. [8] shows a subset of the generated particle spectrum for the parameter point 
PP1. The variation of the particle masses is drawn with a white box in case of 
trivial matching and with a black box in case of implemented threshold corrections. 
The gluino, the lightest chargino and the lightest neutralino show the biggest de- 
pendency with 65 %, 83 % and 85 %, respectively. With the implemented threshold 
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Fig. 7: Dependency of the gluino mass at <5fi x = 3 TeV on the matching 
scale Te 6 ssm for parameter point PP1. The circles show the behavior without 
threshold corrections and the squares with corrections for gi and yf. 



corrections, one finds a reduction of the variation down to 0.1-4 %. The biggest 
decrease is found for the gluino, the lightest chargino and the lightest neutralino 
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whose remaining variation is of the order of 0.5 %. This is because the gluino mass 
is very sensitive to g^, as can be seen from the RGEs of the soft parameter M 3 
[8]. Note that the remaining variation of the particle masses is due to 2-loop and 
iterative effects, which are of the order of up to 4 %. 
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Fig. 8: Particle spectra for parameter point PP1. The white and the black 
boxes show the variation of the particle masses when the matching scale 7e 6 ssm 
is varied in the interval [\Tq, 2Tq], where To = 1.9 TeV is the geometric average 
of all shown particle masses. The black boxes show the error with threshold 
corrections and the white boxes without. 



Note that the variation of the matching scale is not always a good estimation of 
the theoretical uncertainty. For example the error band of without threshold 
corrections in Fig. [8] is shifted when adding threshold corrections but is not reduced 
in size. As can be seen in Fig. 0the reason for the shift is the nondinear behavior of 
rrif) when T Ef . SSM is varied. In the interval [\Tq, 2T ] the exotic mass without 
threshold corrections happens to have a minimum which leads to an abnormally nar- 
row error band. In contrast mjj 2i shows the typical, approximately linear behavior 
in [|T , 2T ] which leads to much better uncertainty estimation. 
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Fig. 9: Dependency of the first generation scalar exotic masses at Qr x = 
3 TeV on the matching scale Te 6 ssm for parameter point PP1. The circles show 
the behavior without threshold corrections and the squares with corrections 
for gi and yj. In the interval [^Tq, 2Tq\ the exotic mass J^/j without threshold 
corrections has a minimum which leads to poor error estimation, whereas wi^j 
is approximately linear within this interval. 



5.3 Exploration of the E 6 SSM parameter space 
5.3.1 Parameter space 

We now turn from the specific benchmark point PP1 to a fuller investigation of the 
EgSSM parameter space. Since our main focus is the impact of threshold corrections 
we still restrict ourselves to a two-dimensional slice of the 12-dimensional parameter 
space. We choose to keep the Yukawa couplings of the SU(5) 5-plets containing Hu, 
H211 Di and Di the same and generation- independent, thus leaving a single unified 
exotic Yukawa coupling, defined at the GUT scale Mx, 

Ai = A 2 = A 3 = ki = k 2 = k 3 (67) 

between those states and the third generation singlet. In addition we fix 

s = jj! = m H > = mg, = 10 TeV, Bjj! = 0. (68) 

The remaining input parameters are tan/3 and A3, and we scan over them in the 
range 

tan p e [2,45], A 3 e[0,3]. (69) 

Note that we are defining A3 at the GUT scale, and in this respect the numerical 
value should not be compared with A in the NMSSM which is often defined near the 
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electroweak scale (EW). By defining A3 at the GUT scale we ensure automatically 
that it is perturbative at all scales between the electroweak and the GUT scaled 

For each choice of (tan (3, A3), requiring electroweak symmetry breaking simul- 
taneously with fl68|) leads to up to four solutions for (mo, M1/2), which we number 
consecutively. In the following (m , M!/ 2 )-plots we will show all of them, where 
solution 1 is preferred over solution 2 in the overlap region. In Figs. [TSHTIJl we will 
show only the first solution for simplicity. 

Since our input parameters (tan /3, A3) have less direct physical meaning than 
the output parameters (m , M 1( / 2 ), we display many results in the (m , Mi/ 2 ) plane. 
To better understand the connection between these two sets we show in Fig. [Till the 
mapping of (tan 0, A3) — > (mo, M1/2). In the top panel we show as a color contour 
plot how tan/3 varies over the (mo,Mi/ 2 ) plane, while on the lower plot we do the 
same for A 3 . Fig. [10] thus allows to read off the original input values for (tan/3, A 3 ) 
for any given point in the subsequent plots. 

The butterfly shape of Fig. [10] is due to the superposition of the four solutions. 
Most of the (mo, Mi/ 2 ) space is covered by relatively small A3 < 0.4, while larger 
values of lambda are concentrated in a narrow region of the parameter space due to 
the renormalisation group flow leading to a focusing on a significantly smaller range 
of A 3 at the electroweak scale. The tan (3 dependency depends on the solution, but 
in much of the parameter space large values of tan (3 are associated with large mo 
and small Mi/ 2 . 

Fig. [TD]is not exactly symmetric in ±Mi/ 2 because A3 and therefore the effective 
/^-parameter is required to be positive in our analysis. 

5.3.2 Experimental constraints 

At this point we briefly summarize the various experimental constraints on the 
model. A more detailed discussion can be found in [TO] . 

Since the publication of the first studies of the CE 6 SSM [S][H] the Atlas and CMS 
experiments have already placed strong constraints on supersymmetry from the first 
« 5 fb _1 of data gathered from the LHC running at y/s = 7 TeV. With 4.7 fb _1 
of data Atlas have performed searches for squarks and gluinos, by looking for jets 
plus missing transverse momentum and possibly one isolated lepton [431 [44] [45] and 

6 In the NMSSM it is well known that to ensure perturbativity A(EW) < 0.7, while an analysis 
of this in the E 6 SSM [2J found that A 3 (EW) < 0.85 is required for perurbativity up to the GUT 
scale, though this limit depends on the values of the other exotic Yukawa couplings and tan /3. 
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Fig. 10: Illustration of the mapping (tan /3, A3) — > (mo,M^2). The top panel 
shows the input parameter tan /3 varying across the (mo, M^^) plane (of output 
parameters). The bottom panel shows the same for the input parameter A3. 

presented the constraints in a CMSSM interpretationE CMS collaboration have 
also placed similar constraints using the razor analysis [46] with 4.4 fb _1 of data 
[47] and another analysis with the transverse mass variable [45] using 4.73 fb" 1 of 
data. Very recently the ATLAS released updated results with 5.8 fb" 1 data from 
the y/s = 8 TeV run [49]. 

7 Thcy also presented constraints in a simplified model interpretation, however this is less rel- 
evant for our present purposes since the physical spectrum of our model, especially if the colored 
exotics are heavy is far closer to that of the CMSSM than the simplified models. 
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It was estimated in [50j EE] that the effect of these constraints on the CE 6 SSM is 
to place a limit on the gluino mass, close to 850 GeV. Since the squarks must always 
be heavier than the gluino, due to the RG flowjfl squark limits are automatically 
satisfied. However the 8 TeV data has increased the limits on the gluino, in the 
heavy squarks region the limit varies between « 900-1000 GeV. 

However it is also of interest to point out that in certain E 6 SSM scenarios (where 
the bino-like neutralino is not stable, leading to longer cascade decays and an alter- 
native LSP) the gluino limits from these searches may not apply [18] resulting in 
different limits. 

Furthermore, there are limits on the Higgs mass [5TI 152] from analyses using 
between 4.6 and 4.9 fb _1 of LHC data. At 95 % confidence level Atlas exclude a 
standard model Higgs in the ranges 110.0-117.5 GeV, 118.5-122.5 GeV and 129- 
539 GeV, while at the same confidence level CMS exclude 127.5-600 GeV. This 
leaves a narrow allowed window of 117.5-118.5 GeV and a wider one of 123.5- 
127.5 GeV. 

As this paper was being finalized the discovery of a new particle consistent with 
a Higgs boson with a mass of around 125-126 GeV was announced by ATLAS [28] 
and CMS [2S]. CMS report a 5.0a excess at 125.5 GeV and quote a best fit mass 
of 125.3 ± 0.4(stat.) ± 0.5(syst.) GeV while ATLAS find a 6.0a excess at 126.5, 
and using the two channels with the highest mass resolution estimate the mass as 
126.0 ± 0.4(stat.) ± 0.4(syst.) GeV. 

The strongest current constraint on the Z' N comes from a CMS analysis using 
« 5 fb" 1 of data [53], announced as this paper is being finalized and setting a lower 
limit of 2.08 TeV on the mass of the Z' N boson. However one should note that the 
lower limit assumes there are now light exotics available for the Z' to decay into. 
If such are present then the branching ratio to leptons is diluted and the limits are 
weaker. 

Finally the model also contains exotic (colored and inert) fermions and sfermions. 
The inert states are weakly produced and we do not anticipate very large new limits 
to be set by the early data from the LHC, but the colored states should be produced 
abundantly and could be tightly constrained. However the existing analyses do 
not apply to the fermionic diquarks or leptoquarks described in this model since 
they carry half integer spin and are odd under their respective discrete symmetries, 
decaying with missing energy. Thus it is currently unknown what limit the LHC 

8 The soft scalar masses get large contributions (relative to the gauginos) from M1/2, due to the 
renormalization group evolution, which is very different from the MSSM case. 
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has placed on their mass. Here we will only consider cases where these exotics are 
heavy anyway but when discussing existing benchmarks from the literature we will 
not apply any new constraints on them. 

The EgSSM also contains a number of dark matter candidates. It is possible that 
the correct relic density can be achieved entirely from the inert sector neutralinos 
(admixtures of inert singlinos and inert higgsinos) [Mj however such a scenario is now 
in conflict with XENON-100 limits [55] though if one allows for a relic density which 
is too small to explain observations (and therefore would require some additional 
contributions to dark matter) consistency could be achieved. Alternatively another 
possible scenario described in Ref. [5B| is that there is a large splitting between the 
inert states with one at a few keV and another a few GeV, thus giving warm dark 
matter scenarios, which match observation. For cases studied the couplings which 
are assumed to be negligible here are too large to avoid perturbing the RG evolution 
and may be too large for consistency with the constrained model generally. However 
such scenarios have only recently been proposed and the full parameters space far 
from explored so it is simply unknown whether or not this could be applied with 
couplings which would be consistent with the scenarios explored here. 

A third possibility which is known to be consistent with the constrained version 
of the model decouples the singlinos from the rest of the spectrum, giving a new 
contribution to the effective number of neutrinos (which is consistent with observa- 
tion) [57] . The bino is then the dark matter candidate, but the scenario is still very 
distinct from typical models with bino dominated dark matter because the inert 
Higgsinos play a crucial role in enabling the correct relic density to be achieved so 
long as the mass of the lightest inert neutralino is related to that of the lightest 
neutralino by the condition [57] . 

fi Sl « \m^o\ + 10 GeV. (70) 

In the scans we perform here we do not explicitly satisfy this constraint (and instead 
our inert Higgsinos are typically far heavier) as we assume a universality amongst 
all exotic Yukawa couplings. Splitting one generation of all exotics or just inerts 
could solve this problem. Doing so would perturb the renormalization group equa- 
tions and distort the parameter space somewhat, though the qualitative results and 
understanding of the impact that these thresholds can have would not be altered. 
Alternatively one can keep in mind the "keVins and GeVins" scenarios though it is 
currently unknown if this could result in the correct relic density or not. 
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5.3.3 Change in masses from threshold corrections 

Now we look at the significance of the threshold effects throughout the parameter 
space described in Sec. 15.3.11 Like for the point PP1, the largest corrections are 
observed in the gaugino sector for states whose masses are set by Mu2- Fig. [TT] shows 
the shifts in mass for the gluino and lightest neutralino in the (m , M 1 / 2 )-plane|^| We 
compare between results with trivial matching conditions to those with full threshold 
corrections, where in both cases the matching scale used is the "optimal" choice, 
where Te 6 ssm is set to the geometric average of the particle masses. 

We show the absolute change in masses rather than relative change because, 
when the masses are light due to a cancellation, the relative corrections can be very 
large, distorting the plot. However it is important to note that the gaugino masses 
though not shown can be estimated since the low-energy soft gaugino masses Mi 
(which give the dominant contribution to nig and m^o) are proportional M1/2, with 
nig ~ 0.85MX/2, m x° 1 ~ 0.15Mi/2- Fig. [T21 shows the actual gluino and neutralino 
masses after including threshold corrections for comparison. 

Once again we see substantial changes in both the gluino and neutralino mass 
across a wide range of the parameter space with corrections being small only in 
a very narrow region of the space. Around the LHC limit on gluinos we can see 
there can be large corrections of the order of several 100 GeV for the gluino mass, 
therefore these corrections can make a significant impact in determining whether or 
not a point is ruled out by LHC search constraints. The corrections to the gluino 
mass can be positive or negative, depending on whether the masses appearing in the 
logarithmic corrections are larger or smaller than the threshold scaled 

The corrections to the neutralino are also substantial, varying from just below 
zero to just over 200 GeV. They are distinctly correlated with the corrections to the 
gluino, showing close to identical variation, which is because the dominant effect is 
a shift in M1/2. 

With the recent discovery of a Higgs-like boson [28| 129] with a mass ~ 125-126 
GeV and strong constraints on a SM-like Higgs away from the mass of this new state, 
the lightest Higgs mass is a crucial observable in constraining the parameter space 
of SUSY models. The corrections to the lightest Higgs mass from threshold effects 
are not expected to be very large since it's mass is not set by soft mass parameters 

9 Notc that these are the too and -^1/2 values obtained from the spectrum generator when thresh- 
olds are included. Since m and Mi/ 2 are outputs they are affected by the threshold corrections 

and so do not match the values for the point without threshold corrections. 

10 In this case the dominant corrections arise from corrections to #3, see Eq. (|55]h 
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Fig. 11: Change in the gluino and in the lightest neutralino mass when adding 
threshold corrections for s = // = m#> = mgi = 10 TeV, tan/3 £ [2,45] and 
Aj(Mx) = Ki(Mx) € [0,3]. 

at tree level. However even small corrections to the Higgs mass can significantly 
shift the mo and M1/2 values which a particular Higgs mass is compatible with and 
this can dramatically alter how the Higgs mass measurement combines with the 
constraints coming from squark and gluino searches. 

In Fig. [13] the corrections to the Higgs mass are shown across the (mo, M1/2)- 
plane. The corrections vary between around —0.7 GeV to +1.7 GeV. One can also 
see in Fig. [TUthe Higgs mass itself and note that for the parameter set studied here 
the Higgs mass is consistent with the new particle discovered by ATLAS and CMS 
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Fig. 12: Gluino and in the lightest neutralino mass as a function of Mi/ 2 
and mo, where it was varied tan/3 6 [3,50], Aj = k« £ [0,3], s = 10 TeV, 
// = niH' = niff, = 10 TeV. 



in the large mo and M\/2 regions of the parameter space where the corrections from 
the thresholds are largest. 

5.3.4 Allowed parameter space 

With the more reliable predictions at hand, we turn to understanding the overall 
allowed region of parameter space. Dominant experimental constraints come from 
the gluino and the Higgs, see Figs. [12] and [TH We can see that for our choice of s the 
gluino is comfortably above the LHC limit in a large volume the parameter space. 
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Fig. 13: Change in the lightest CP even Higgs mass when adding threshold 
corrections for s = // = m#/ = m H' = ^ TeV, tan/3 E [2,45] and Xi(Mx) = 
Ki(Mx) E [0, 3]. We show only points where the change is less than 5 GeV. 
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Fig. 14: Lightest CP even Higgs mass as a function of M]/ 2 and mo, 
where it was varied tan/3 £ [3,50], Aj = Kj £ [0,3], s = 10 TeV, // = = 
m^, = 10 TeV. 

In the Higgs contour plot Fig. [14] we also see a very large volume of the parameter 
space where the light Higgs mass is within the narrow window on it's allowed mass 
from the LHC discovery. Nevertheless in fact the combination of the two constraints 
creates very stringent limits of this slice of the parameter space. 
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To show the impact more precisely we now plot in Fig. [15] the valid and invalid 
parameter space including all relevant constraints, both in terms of (tan/3, A3) and 
of (m,Q, M1/2). While it is beyond the scope of the current project to take account 
of full experimental likelihoods or include precise 95% confidence limit contours, to 
provide a guide as to the parameter space which evades LHC constraints we define 
valid points as those satisfying, 



which are based on the discussion of Sec. 15.3.21 

One finds that the lower bound on m g and the bounds on are most constrain- 
ing. This is because the gluino is driven lighter than the sfermions by the renor- 
malization group flow with, M 3 ~ 0.7Mi/ 2 ~ 200-1400 GeV, while for the Higgs 
we now have a very narrow range of allowed masses as a result of Higgs searches 
that ultimately lead to the recent discovery. Here we find the lightest Higgs mass 
is around m/ tl ~ 118-128 GeV [8], as shown in Fig. [T4] including a substantial re- 
gion of parameter space where the CEgSSM predicts a lightest Higgs with a mass of 
123.5-127.5 GeV, consistent with the discovery. 

While the allowed region appears strongly restricted in (tan/3, As)-space, it ap- 
pears reasonably large in the (mo, M\ 12) -plane. The reason is the non-trivial map- 
ping between these two spaces, as illustrated by Fig. [TO]. The valid green region in 
the (tan /3, As)-space of Fig. [15] is mapped onto different wings of the butterfly in the 
(m , Mi /2 )-plane. 

Fig. [T5] also shows that even without imposing experimental constraints there 
are upper and lower limits on all the parameters. Outside the colored regions no 
simultaneous solution to electroweak symmetry breaking and unification without 
tachyonic masses can be found. In particular, there is an upper limit on tan/3 
at around 45. Here the EWSB minima become unstable giving a tachyonic Higgs 
mass, as shown in Fig. [16] for both the uncorrected Higgs mass and for the threshold 
corrected value, which converge in this region. At such large tan /3 the bottom quark 
Yukawa coupling is much closer to the top quark Yukawa coupling, and in such a 




123.5 GeV < m hl < 127.5 GeV, 
m g > 1 TeV, 
m h > 300 GeV, 
m^o,m-± > 200 GeV, 
m z , > 2.1 TeV, 



(71) 
(72) 
(73) 
(74) 
(75) 
(76) 
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Fig. 15: Exclusion plot of the EgSSM parameter space. 



situation generating a large splitting between the Higgs masses required for correct 
EWSB is difficult. A simple tree-level condition for keeping m\ > at large tan/3 
can be derived in the MSSM [58] and in the E 6 SSM, 

MSSM: m 2 Hd - m 2 Hu > m|, (77) 
E 6 SSM: m\ A - m 2 Hu > m% + 9 -f^ Q (N H2 - N Hl ). (78) 

The E 6 SSM condition is stricter since it contains an additional positive term involv- 
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Fig. 16: Second lightest CP even Higgs mass as a function of tan f3 for s 
fjf = m H > = m s > = 10 TeV and \i(M x ) = Ki{M x ) = 0.2 
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(79) 



5.4 Dependency on the survival Higgs parameters 

The inclusion of threshold corrections not only reduces the unphysical matching- 
scale dependence. It also leads to a physical dependence on the survival Higgs 
sector. Without thresholds the spectrum does not explicitly depend on the survival 
Higgs masses as the renormalization group equations decouple. With threshold 
corrections the survival Higgs parameters appear in the threshold corrections to the 
gauge couplings, and thus affect both gauge coupling unification and the low-energy 
mass spectrum. We will now study these two aspects in a simplified setting by 
taking a single scale 



mg,. 



50) 



Figs. [T7] and [TH] focus on gauge-coupling unification. We recall that in our ap- 
proach Mx is defined by the intersection of g\ and g-i , while #3 is determined by its 
running from its low-energy measured value. Fig. [T7] shows that for the parameter 
choice PP1, exact unification can be achieved by adjusting m surv = 200 GeV, while 
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The term is positive so long as s > v, which is always the case due to limits on the Z' mass. 
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Fig. 17: Splitting between g\ and 53 at Mj for two different choices of m s 
for PP1. 



Of course in a full GUT model there will be GUT scale threshold corrections. 
Nevertheless Fig.[T7]makes clear that requiring a precise splitting between the gauge 
couplings from any given GUT threshold correction will fix the survival Higgs scale. 

Fig. [18] shows gauge coupling (non-)unification for the entire parameter space 
defined by Eqs. (I67j) - (l69|) except that tan/3 = 10 is fixed and m surv is varied in the 
range 

WW G [0.1, 1000] TeV. (81) 

Like in Fig. [TTJ we find that the deviation (g 1 — #3) at M x depends substantially on 
the survival Higgs masses, and for any given lambda there is a survival Higgs scale 
which allows exact unification or any splitting required by GUT thresholds between 
and O(0.1). 

Figs, dn] and [2U] focus on the influence of m surv on the mass spectrum. In Fig. [TH] 
we plot the variation of the gluino mass across the (m surv , A3)-plane and see that the 
dependence on the survival Higgs masses is extremely weak. The figure also shows 
that the lightest neutralino mass is affected slightly more, in particular in regions 
where it is rather heavy (~ 250-300 GeV). 

Typically the survival Higgses and Higgsinos are too heavy to be observed at 
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Fig. 18: Splitting between gi and at Mj in the (A3,m surv ) parameter 
space 



the LHC, since they are only weakly produced and have the couplings of a 4th 
generation lepton with current mass limit on the mass of 100.8 GeV. 

For cases where they are not observable and where we do not constrain the 
splitting between the gauge couplings at the GUT scale, we can therefore consider 
them as an additional error in the theoretical calculation. This is the approach we 
now take here. To look at this in detail we again turn to our test case PP1. 

Fig. [2Q] shows the dependency of the particle spectrum on the new model pa- 
rameter m surv for PP1. It is instructive to compare the result to Fig. [Sj where the 
sensitivity to the matching scale Te 6 ssm with and without threshold corrections is 
shown. In the comparison one has to bear in mind that the m surv -dependence is 
physical and that we vary m surv in a much larger interval than Te 6 ssm- 

One finds a significant m surv - dependence only for masses which were also strongly 
sensitive to the choice of matching scale without threshold corrections. This is be- 
cause a variation of Te 6 ssm changes the size of all the logarithmic contributions to the 
threshold corrections, while the variation of m surv changes a subset of contributions 
to the thresholds. However, we also find that after including threshold corrections 
the m surv -dependence is generally larger than the remaining matching- scale depen- 
dency. 
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Fig. 19: Gluino mass m„ and mass of the lightest neutralino m^o in the 
(A3,m surv ) parameter space 



5.5 Benchmarks in the literature 

As a final application of the threshold corrections presented here we now study the 
impact they have on previously published benchmarks in the model which have 
appeared in [HI El CD] • The 1 TeV limit for the gluino which we require to be consis- 
tent with Atlas and CMS searches is above the masses of the all light benchmarks 
[H [9] , while heavy benchmarks [10] were chosen to be safe from these limits. How- 
ever since the threshold effects presented here have a very large impact on the gluino 
mass, clearly they play an important role in determining whether or not a particular 
CE 6 SSM point is excluded. We will therefore study what impact the threshold ef- 
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Fig. 20: Dependence of the particle spectrum on the survival Higgs scale 
TTisurv = 0.1-10 TeV for parameter point PP1. 



fects have on all these points and in particular whether it changes their experimental 
statusH 

The benchmarks in Refs. [H |9] were selected to have very light gauginos or be 
close to the lower limit on m (for that choice of singlet VEV s), where one can get 
lighter exotic sfermions. Below this lower limit on m the inert Higgs scalars become 
tachyonic due to the large U(1)n .D-terms which give a negative contribution to the 
mass. 

As a result of this the threshold corrections can push some points into the region 
with tachyonic masses. Since we are applying the threshold corrections iteratively 
this creates a problem. It may be that such points in fact do not contain tachyonic 
solutions but merely the first step in the iteration jumps to that region and sub- 
sequent steps, if they could be performed, would lead to a self consistent tachyon 
free solution. In this case we try pushing the point into a convergent iteration by 
adjusting the thresholds corrections in the early steps. If this does not work we then 
vary the survival Higgs mass to see if this can lead to tachyon free solutions. 

The impact on the gluino mass of these threshold corrections, for all the points 

12 All but one of the light benchmarks have a Z' mass substantially below the limit discussed 
in Sec. 15.3.21 Light exotics, like singlinos, which are necessarily present in the EgSSM, will give a 
substantial reduction in the limit. This is probably not sufficient to evade the limits, but nonethe- 
less testing against constraints from gluino searches provides another important limitation on the 
viability of these scenarios. 
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nig / GeV 


BMA 


BMB 


BMC 


BMD 


BME 




without thresholds 


336 


330 


353 


327 


338 




including thresholds 


224 


269 b 


260 


230 


203 b 




nig / GeV 


BM1 


BM2 


BM3 


BM4 


BM5 


BM6 


without thresholds 


350 


673 


362 


642 


338 


805 


including thresholds 


322 


613 a 


275 


423 a 


190 


825 


m g / GeV 


HBM1 


HBM2 


HBM3 


HBM4 


HBM5 




without thresholds 


984 


1352 


1659 


1129 


1001 




including thresholds 


1090 


1494 


1886 


827 


1067 





Tab. 4: Comparison of the originally reported gluino mass and the gluino mass 
including threshold corrections for previously published EgSSM benchmarks, 
where we have chosen // = tuh' = m B> = 10 TeV for all points except those 
marked a where we used a 1 = ran 1 = fngi = 100 TeV to evade tacbyonic 
problems with convergence and b where they were increased to 10 4 TeV. 



is shown in Tab. H] for y! = mH 1 = mgv = 10 TeV, unless marked with a or b where 
the survival Higgs masses were varied until a self consistent solution could be found. 

But while the survival Higgs masses can have important effects, all light bench- 
marks are left with a gluino mass substantially below the limit. Therefore we can 
confirm that all of these points have been excluded by the LHC even when threshold 
effects are taken account of and the survival Higgs masses are varied without regard 
to maintaining gauge coupling unification. 

The situation is more optimistic for the new heavy benchmarks proposed in [10] . 
These benchmarks were chosen to satisfy the latest LHC limits and also have the 
right relic density. For the benchmarks HBM1-HBM3 there are significant changes 
to the gluino mass as shown in the third column of Tab. H] pushing them further 
away from current constraints. However beyond this all qualitative features remain 
unchanged with the Higgs mass staying in the Higgs signal range, though it is 
modified to 125 GeV for HBM1 and HBM3. The light inert Higgsinos, fixed to 
be light to achieve the correct relic abundance density, remain light and can be 
tuned carefully to match the relic abundance density measurement exactly without 
perturbing the rest of the spectrum. HBM2 and HBM3 also have light exotic quarks 
and these also stay light, with their masses reduced by approximately 4 % and 6 % 
respectively. Therefore the essential features of these benchmarks are not changed 
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by the inclusion of threshold effects. 

However for HBM4 and HBM5 the situation is more dramatic. These points 
have a very large singlet VEV, 50 TeV and 100 TeV respectively and in such cases 
while light gauginos can always be obtained, this is only achievable for a very narrow 
range of the input exotic Yukawa couplings. Consequentially when one varies these 
Yukawa couplings the masses can change a lot and the modification in these values 
(at low energies) from threshold corrections can result in very large changes to the 
soft masses. For HBM4 although a gluino mass of the same order is maintained the 
sign of M1/2 is actually changed, while for HBM5 a positive Mi/ 2 is obtained but 
the stability of the result is not clear as scale variations lead to huge changes in the 
mass. 

However these large changes in the physical spectrum are really an artifact of 
the choices for input/output parameters. One can instead try to match the values 
of mo, M\j2 and then observe changes in the exotic Yukawa couplings. This would 
typically lead to a similar spectrum, though A3 affects the light Higgs mass and 
therefore the light Higgs mass is shifted. 

6 Summary and conclusions 

In this paper we have improved the prediction of low-scale quantities from high-scale 
parameters in the Exceptional Supersymmetric Standard Model (E 6 SSM). Full spar- 
ticle threshold corrections to the gauge couplings have been calculated by match- 
ing the E 6 SSM to the SM. Similarly, the low-scale DR E 6 SSM Yukawa couplings 
have been calculated at the loop level directly from mass measurements of the SM 
fermions and raw Full expressions for these corrections are provided, and we have 
implemented them into an improved spectrum generator for the constrained EgSSM. 

Using the spectrum generator, we have studied the impact of these corrections 
in detail. Both for the test point PP1 and in a parameter scan we found a dra- 
matic reduction in the scale dependency of many masses, in particular of the gluino 
mass which is important for setting limits on the model from collider constraints. 
Also charginos, neutralinos, squarks and exotics receive large corrections. Even the 
smaller threshold effects for the light Higgs boson mass have a big impact on the 
allowed regions of parameter space due to our knowledge of the allowed values of 
the light Higgs mass, given recent LHC results on Higgs searches and the discovery 
of a new particle. 

The new spectrum generator with the implemented corrections allows us to draw 
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firmer conclusions about the allowed parameter space of the model. In line with 
previous results, we found a substantial region of parameter space where the mass of 
lightest Higgs is compatible with the new boson discovered at the LHC and the limits 
from collider searches for squark, gluinos and Z' bosons, are all satisfied. However 
the threshold corrections imply a significant change in the high-scale parameters 
where this is achieved and can alter the way these constraints combine. 

An interesting consequence of the threshold corrections is a dependency of the 
mass spectrum on the survival Higgs fields, which are included in the model to assist 
gauge coupling unification. In previous studies these states simply decoupled from 
the rest of the spectrum. With the full sparticle threshold corrections included, we 
found that in most of the spectrum the dependency of the spectrum is rather weak. 
However if the dependence on the unknown survival Higgs masses is viewed as a 
theoretical uncertainty, this uncertainty is larger than that from the remaining scale 
variation, though much smaller than the uncertainty of earlier studies neglecting 
threshold corrections. 

If the survival Higgs bosons are taken seriously as physical fields, the threshold 
corrections allow to fix their masses by the requirement of gauge coupling unification. 
Indeed we found that the survival Higgs masses can always be chosen such that 
gauge coupling unification is valid, up to hypothetical GUT threshold corrections, 
with required survival Higgs masses in the multi-TeV region. 

Finally we looked at how the threshold corrections affect benchmark points that 
have previously appeared in the literature. Many benchmarks proposed in earlier 
papers appeared to be ruled out by gluino searches. We confirmed all of these to 
be ruled out once the threshold corrections are included. The benchmarks pro- 
posed more recently to be consistent with experiment receive significant corrections. 
Nonetheless these benchmarks are still experimentally viable after including the 
threshold corrections. For two benchmarks, HBM4 and HBM5, the stability of the 
predictions is poor due to fine tuning of input parameters required to get the very 
hierarchical spectrum these benchmarks illustrate. 

In conclusion, our study shows that it is possible and valuable to take into account 
higher-order corrections to the high-scale-low-scale connection even in complicated, 
non-minimal supersymmetric models such as the EgSSM. The qualitative features 
of the CEgSSM itself have not been changed by the threshold corrections. Still, the 
model is theoretically attractive, predictive and viable. As more LHC data comes 
in, distinguishing between different supersymmetric models and hypotheses of GUT- 
scale physics becomes more relevant, and precise spectrum generators like the one 
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presented here will be helpful. 
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A EgSSM covariant derivative and GUT relations 



The Eq covariant derivative reads 

D^d^ + igoAlf" (a = l,...,78), (A.l) 

with a single gauge coupling go, where T a are the generators and A ^ are gauge 
fields in the adjoint representation (78) of the E§. This covariant derivative can be 
decomposed in terms of the SU(3) C , SU(2) L , U(l)y and U(1)n sub-groups 

t - Y N 

D, = d, + ig 3 T a G; + ig 2 - ■ W, + \gy-B, + ig N -Z^ + ■■■ (A.2) 



All generators in (1A.2|) are normalized such that the quantum numbers are the same 
as in Tab. [TJ This, together with the condition Tr T a T a = 5 ab /2, yields a GUT 
relation between the couplings 

93 = 92 = 9i= g[ =■ go, (A.3) 

where 



9i-=y^9Y, g[:=VA0g N . (A.4) 

At the electro weak scale the generator Q of the unbroken U(l) em is then given by 

«=7 + T (A ' 5) 



B E 6 SSM mass eigenstates 

B.l Higgs sector 

We write the Higgs bosons in the E 6 SSM as 

H P i = \ J™ ) > s i> (B- 1 ) 
\ H piJ 

where % — 1, 2, 3 is the generation index and p = 1, 2 denotes the down (p = 1) and 
up (p = 2) Higgs bosons. When the SU(2) L x U(1) Y x U(1) N is broken to U(l) em , 
the third generation Higgs bosons get a vacuum expectation value 

( H?*\ f^ + ^e H° + i 3m if ° , 
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23 




H. 



23 



23 , 



S3 — > — ^= + 9fte S3 + % Sm 
V2 



Furthermore we define 

tan (3 : = 



V2 



tan i 



2.s 



sin 2(3, 



v/2' 



(B.3) 
(B.4) 

(B.5) 



From the real parts of H® 3 , if 23 and S3 we construct three CP even Higgs bosons. 
The diagonalization of the CP even mass matrix is done in two steps. At first we 
transform 3fteif° 3 , 5teif 23 and JieS^ into intermediate states (hi, h 2 , h 3 ) via 

1 

he 



\ 



+ 



v/2 



5?eS 3 y 

where the mixing matrix C/mssm has the form 

' cos (3 — sin j3 0^ 



''2 

V s / 



(B.6) 



U 



MSSM 



sin 

sin (3 cos (3 
I) 



(B.7) 



In the basis of h = (h 1; h 2 , h 3 ) the Lagrangian for the CP even Higgs masses reads 



£cvcn — — ^ ^MMSSMh, 



(B.8) 



where the matrix Mmssm is non-diagonal in general. Note, that the above trans- 
formation is analogous to the MSSM, where the mixing angle in ?7mssm is (3. In 
a second step we diagonalize the mass matrix Mmssm by the unitary matrix Ue 6 - 
The resulting CP even Higgs mass eigenstates are labeled h = (hi, h 2 , h 3 ) T . The 
diagonalization transformation reads 



h = U Ee h, 



M E , = U E M MSSM Ul 



(B.9) 



where Mg 6 is diagonal. From the gauge eigenstates A = (*<smHl 3 , QmH^, Sm Sj 
we construct three CP odd Higgs boson mass eigenstates A = (Ai,A 2 ,A 3 ) T = 
(G°, G', A°) T via 



V2 



A := U A A, 



M A = U* A MU\, 



(B.10) 
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(B.ll) 



where Ma is diagonal. The mixing matrix Ua is parametrized as 

^ cos (3 — sin (3 ^ 

Ua = — sin /3 sin — cos ,5 sin cos0 
i sin j3 cos cos /3 cos sin 0y 

The charged Higgs and goldstone bosons (H i ) = (G^ H ± ) T are constructed from 

■± 

cos /3 — sin (3 
sin /3 cos j3 



the gauge eigenstates if^ via 

_1_ _L -L- x _L 

Hf = Ut 3 H% (z,j = l,2), where f/± = ( ^ (B.12) 



B.2 Inert Higgs sector 

The first two generations (z = 1, 2; p = 1, 2) Higgs doublets in Eq. ( IB.lj) are called 
inert Higgs bosons. For each generation i = 1, 2 we mix the fields H^, H®* to mass 
eigenstates h° ik with an unitary matrix U^ crt via 

n , n • s f H?f\ n - / cos#,° sin# 7 °\ 



Here k, I = 1,2 enumerates the mass eigenstates and we neglect inter-generation 



mixing. Furthermore for each generation i = 1,2 we mix the fields H u , H^* to 



mass eigenstates h ik with an unitary matrix U^ ert via 



K k = ( ) , ^ t = I ~t:+ "";v i • (B.i4) 



H 2i 




Here /c, / = 1,2 enumerates the mass eigenstates and we neglect inter-generation 
mixing. Furthermore for each generation i = 1,2 we mix the fields H® iL , H^n to 
mass eigenstates ip^ with an unitary matrix Z via 



< = Z ln 2 a > Z = 4 L ■ (B-15) 





The Majorana mass eigenstates are then defined as 



h^=[3r). (B.16) 



For each generation i — 1,2 we combine the fields i? Xii/ , to mass eigenstates 

h~ via 

V = I • (B-17) 

\#2iL 

50 



B.3 Survival Higgs sector 

We mix the neutral survival Higgs bosons H'°, H' * to mass eigenstates h'° k with an 
unitary matrix t/ s ° urv via 

/0 ( H'°\ / cosO' sm6'°\ 

h k = (U surv ) kl ^j, U sarv =^ ffQ cos9 , j. (B.18) 

We mix the charged survival Higgs bosons H' , H' to mass eigenstates h k with 
an unitary matrix U^ TV via 

, + s f H'-\ . ( cosfl'* sinfl'A 

^ = (^4H; ^ cos^j' (B19) 

The survival higgsinos obey the same mixing as the neutral and charged inert hig- 
gsinos in Sec. IB. 21 We write the mass eigenstates as h !± and hf (i = 1,2). 



C E 6 SSM self-energies 



C.l W boson 

The W ± boson 1PI correlation function is decomposed into a transverse and longi- 
tudinal part as follows 



p2 



TlwwAP 2 ) - ^-^-ww,l{p 2 



(c.i; 



where the transverse part is in the E 6 SSM given by 

3 2 



92 



i=l k=l 



+ \A Q {m H± ) + U (m G ±) + m 2 w £ {U^*) 2 B (p 2 



i=i 



2 uv 2 

3 2 

EE \ U * U n + Uf 2 U%\ 2 B 22 {p\m Av m Hf ) 
i=i j=i 

3 2 

EE 

i=i j=i 

2 



m hi ,m w> 



E^ 6 (^ SSMf/ fi-^ SSMf/ 3 



/c=i 



B 22 (p 2 ,m Ai ,m H± ) 



E{^(I^| 2 + I^| 2 ) 
U=i ^ 
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- 2Z ll Z l2 rni+mto B (p , m~ h+} m~ h o ) 

i li i li 

2 

- 2 { cos2 (^° + e t) (b 22 (p 2 , m^m^) + B 22 (p 2 , m h ± , m^)) 

+ sin 2 (fl° + Of) (b 22 {p 2 , m h ±,m h o 2 ) + B 22 (p 2 , m h ±, m^)) } 

+ EU(I^i| 2 +I^ 2 | 2 ) H^m-^m^) 
1=1 ^ 

- 2Z' ll Z' a mf lf+ mf l pB {p l , m~ h , + ,m~ h , ) 

2 

- 2 ^ { cos 2 (#'° + 0'*) (s 22 (p 2 , 7tv±, m h ,o) + 5 22 (p 2 , m h ,±, m h ,o)' 

i=l 

+ sin 2 (# /0 + 7± ) (§ 22 (p 2 ,m h ,±,m h ,o) + B 22 (p 2 ,m h ,±,m hf o)^ } 

i /j 2 9n ( { Nn ia . N H \ 2 . 2 2 2 
+ 477v-|-<! I — 1 — J sin /3cos /3-B (p ,m H ±,m Z ') 



N H 23 -2 o N H 13 „ 9 



sin /S — ^ cos /3 B G {p ,m G ±,m z 



2 ' 2 

- (85 22 (p 2 , m w , 0) + 4p 2 B (p 2 , m w , 0) 

- {(4p 2 + m| + m^c 2 , - m|4} i?o(p 2 ,m z ,m w ) 
-8c 2 w 5 22 (p 2 ,m z ,m w ) 

+ ]T | ^iV//f(p 2 ,m u ,m d ) - ]T 2Nlw 2 ij B 22 {p 2 ,mu. 

fu/f d I ij=l 

i 6 2 r 

+ ~ ^2 ^2 1 fijwH{p 2 , m^o,m^+) + 2g ijW m^om^+B (p 2 , m^o, m^+) 

(C.2) 

Analogous to [5H] the summation Ylf u /f d * s over Q uar k and lepton doublets and 

(w fij ) = ( CuCd CuSd ) . (C.3) 

The neutralino-chargino-W-boson couplings are given by 

fijw = \ a x°xlw\ 2 + \hlx+w\ 2 i 9i 3 w = 2 Ke (b*^+ w a m + w ) , (C.4) 
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where the Feynman rule for the neutralino-chargino-W^ vertex is written as —i r ) ll {aVL+ 
bVn). The nonzero couplings in the E 6 SSM are the same as in the MSSM 



a i>°i>+w 



The couplings to mass eigenstates for an incoming neutralino Xi are 



a x°xjw 



while for an incoming chargino xt the couplings read 



a x°ixt w 



C.2 Fermions 

We decompose the fermion 1PI correlation function as 

r ff -( P ) = f (p L T L fI (p 2 ) + P R rf f -( P 2 )) + P L T l ff -( P 2 ) + P R r r ff -{p 2 ) 

and then define the fermion self-energy to be 



) : = 2 \ m f 
In the E 6 SSM it is given by 



m t 



-y-<j Bxip^mg,™^) + B x {p 2 ,mg,m i2 ) - \h + 31n £| 



sin(2fl t )^ {B (p 2 
m t 



,m~ g ,m 



fj - B Q (p 2 ,mg,m i2 )) 



Vt 



i=l 



\ -t* [ Bl ( p2 ,rn t ,m h J + B (p 2 ,m t ,m hi )] 



(C.5) 



(C.6) 



(C.7) 



(C. 



T ff + r ?/^ 2 )] + V W) + T W) } ■ ( c - 9 ) 



+ B^ [Bi (p 2 ,m t ,m Ai ) - B (p 2 , m t , m A , 
(Vb s l + y1 c l) B i{p 2 ,m b ,m H +) + (g 2 + y^c 2 + y^s^B^p 2 ,m b ,m w ) 



, 2 2 



^(p 2 , m b , m H +) - B (p 2 , m b , m w ) 



(ee 4 ) 2 (5 + 31n^ 



+ 



9l 



Sl + 9l) B.(P 2 , m„ m z ) + 4 ftlSla B„(p 2 , m„ m z ) 



i=l j = l 



fuu B i (P > m x° > m % ) + ft* TT^o (p , , mi. ) 



m t 
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+ £ 



m- 



+ 9% 



2 s 



'ilb. 



m t 



tR 



B 1 (p\m t ,m z ,) + N tL N tR B (p 2 ,m t ,m z , 



(CIO) 



The matrix elements and B_^ are defined as 



.4 



■fi 



B 



fi 



(U M ssu)2k(UEe)tk if / is up-type. 
{Umssm) ik(UE 6 )*k if / is down-type, 

(Ua)* 2 if / is up-type, 
{U a)*i if / is down-type. 



(C.11) 



(C12) 



In analogy to [59] the Feynman rules for the Xiffj couplings are written as —i(dPL + 
bVn) and we have defined 



iffi 



I I 2 _l_ \h \ 2 



g tffj =2Me(b* Mf a Mfj 



(C.13) 



In the gauge eigenstate basis ip Q ,ip + one has 





_ 9Yy 


Hvk - 


9Yy 




(C.14) 


huh 


= V2g 2 ri L , 


a i>?du L 




= 92, 


(C.15) 


a ^dd L 


= H%dd~R = 




u dR 


Vd, 


(C.16) 


a i>°uu L 


= b i>° 4 uu R = 


~ a i>iud L = 


- a i>+du R = 


= Vu, 


(C.17) 


HlffR 


V2 fR1 


Huh = 


y/2 ,Li 




(C.18) 



where the quantum numbers Yf/2, Nf/2 and r/ are listed in the Tab. [TJ The 
couplings to the mass eigenstates xl an d xt are obtained by the rotations 



a x°ff ~ N h a ^fh 
a xtff' = V i*j a i>+ff" 



b x°ff ~ N v b tyff> 



J xjff 



(C.19) 
(C.20) 



To obtain the couplings to the sfermion mass eigenstates one rotates these couplings 
(both a- and 6-type) by the sfermion mixing matrix, 



a xff{ 
l xfh, 



-Sf c f 



l xff' L 
l xff' R < 



(C21) 
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The self-energies for the other fermions are obtained from T, t by the substitutions 



^ 2 )=^ 2 )U 53 =o> 

^(P 2 ) = ^(P 2 )U W . 
In Eqs. (jC.2p and fIC.lOp we use the loop functions as defined in 



(C.22) 
(C.23) 



D EgSSM counterterms 

If not otherwise stated we renormalize the W boson and the SM fermions in the 



on-shell scheme. The corresponding on-shell and DR counterterms are 



5m 



2, on-shell 
W 



Eell 



WW,T\ m W) 



5m 



2,DR 

w 



9fJeIT 



ww,T\jn w j 



5mf- shcl1 



fteZ f (m 2 f ), 



5m™ = sfteE/(m3 



(D.l) 
(D.2) 



where the self-energies 1Iv^w;t and E/ are given in Eqs. (jC.2p . ( IC.lOp . ( 1C.22|) and 
(1C.23p . Using Eqs. (1C2|) and (IC.lOP one can derive the following divergences for the 
5m\y and 5m f counterterms in the E 6 SSM 



5m w 



m w 



A 11 2 1 2 2 

J^)Y 9 * + 2 9y + 29n 



H 13 



cos 2 j3 



H 2S 



2 J V 2 

A 2 - 3y 2 sin 2 - 3y 2 cos 2 - y 2 T cos 2 } , 



sin 2 (3 



(D.3) 



m fu 



Sm f d 



A 



(4tt) 



A 



-2glC 2 (r(f u )) su{3)c 



' 2 2 

l9 2 ~ 9n 



N 



LJu 



-9 2 v 



2 
3 



-2g 3 C 2 (r(f d ))su(3) c - ~ 9 N 



N LJd 



-9 2 y 



RJd 



3 o 2 

+ n y h + n y U 



RJu 



3 1 
+ n y i + o y l } ' 



(DA) 



(D.5) 



Here C2(r)su{N) is a representation invariant of the representation r of SU(N) and 
defined by 

N 2 - 1 



C 2 (iV) 



SU(N) 



2iV 



62(G) SU(N) = TV 



(fundamental representation N). 
(adjoint representation G). 



(D.6) 
(D.7) 
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Especially we have C2(rf u )su{3) c = C < 2( r / d )s , i/(3) c = 4/3. Furthermore it follows from 
the one-loop /3 functions (1D.9P (ID. 12[) that the divergences of Sgi in the SM and the 
E 6 SSM are given by 



(D.8) 

(D.9) 
(D.10) 

(D.ll) 
(D.12) 







A Ft- 
a ft 2 








9i 


A 


" (4tt) 2 2 9i ' 








oEgSSM 
P3 


= 0, 




/3f M = 


-7, 


oEeSSM 


= 4, 




P2 — 


19 


^E 6 SSM 


3 oEeSSM 

- 5 Py ~ 


48 

T' 


/3 1 SM = 


5^ 


o/E 6 SSM 
3 1 


1 oEeSSM 

- 40^ " 


47 
" T' 
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To' 



The divergence of £ tan/3 in the E 6 SSM is 

S tan /? 



tan /? 



A A 



tan / 



(47T) 



from which we obtain the one-loop RGE for tan f3, 

d tan f3 tan f3 



dt 



(4vr) 2 
3 



Aan/3) 
1 



r> ./ ~ 2 , 1 2 ° 2 



iV 



h 23 



(D.13) 



(D.14) 
(D.15) 



needed so that we consistently input tan/3, defined at a fixed scale, when we vary 
the matching scale where the Yukawas are calculated. 



56 



References 

[1] P. Athron, J. P. Hall, R. Howl, S. F. King, D. J. Miller, S. Moretti and R. Nev- 
zorov, Nucl. Phys. Proc. Suppl. 200-202, 120 (2010). 

[2] S. F. King, S. Moretti, R. Nevzorov, Phys. Rev. D 73 (2006) 035009. 

[3] S. F. King, S. Moretti, R. Nevzorov, Phys. Lett. B 634 (2006) 278. 

[4] J. E. Kim and H. P. Nilles, Phys. Lett. B 138 (1984) 150. 

[5] U. Ellwanger, C. Hugonie and A. M. Teixeira, Phys. Rept. 496, 1 (2010) 
[arXiv:0910.1785 [hep-ph]], M. Maniatis, Int. J. Mod. Phys. A 25 (2010) 3505 
[arXiv:0906.0777 [hep-ph]]. 

[6] P.Fayet, Nucl.Phys. B 90 (1975) 104; H.P.Nilles, M.Srednicki, D.Wyler, 
Phys.Lett.B 120 (1983) 346; J.M.Frere, D.R.T.Jones, S.Raby, Nucl.Phys.B 
222 (1983) 11; J.P.Derendinger, C.A.Savoy, Nucl.Phys.B 237 (1984) 307. 
M.I.Vysotsky, K.A.Ter-Martirosian, Sov.Phys.JETP 63 (1986) 489; J.Ellis, 
J.F.Gunion, H.Haber, L.Roszkowski, F.Zwirner, Phys.Rev.D 39 (1989) 844. 

[7] F. del Aguila, G. A. Blair, M. Daniel, G. G. Ross, Nucl. Phys. B 272 (1986) 
413. 

[8] P. Athron, S. F. King, D. J. Miller, S. Moretti and R. Nevzorov, Phys. Rev. D 
80, 035009 (2009) [arXiv:0904.2169 [hep-ph]]. 

[9] P. Athron, S. F. King, D. J. Miller, S. Moretti and R. Nevzorov, Phys. Lett. B 
681, 448 (2009) [arXiv:0901.1192 [hep-ph]]. 

[10] P. Athron, S. F. King, D. J. Miller, S. Moretti and R. Nevzorov, arXiv: 1206.5028 
[hep-ph] . 

[11] R. Howl and S. F. King, Phys. Lett. B 652, 331 (2007) [arXiv:0705.0301 [hep- 
ph]]. 

[12] R. Howl and S. F. King, JHEP 0801, 030 (2008) [arXiv:0708.1451 [hep-ph]]. 

[13] R. Howl and S. F. King, JHEP 0805, 008 (2008) [arXiv:0802.1909 [hep-ph]]. 

[14] J. Erler, P. Langacker, S. Munir and E. Rojas, JHEP 0908, 017 (2009) 
[arXiv:0906.2435 [hep-ph]]. 



57 



[15] R. Howl and S. F. King, Phys. Lett. B 687, 355 (2010) [arXiv:0908.2067 [hep- 
ph]]. 

[16] F. Braam, A. Knochel and J. Reuter, JHEP 1006, 013 (2010) [arXiv: 1001.4074 
[hep-ph]]. 

[17] J. P. Hall, S. F. King, R. Nevzorov, S. Pakvasa and M. Sher, arXiv:1109.4972 
[hep-ph] . 

[18] A. Belyaev, J. P. Hall, S. F. King and P. Svantesson, arXiv: 1203.2495 [hep-ph]. 
[19] T. G. Rizzo, Phys. Rev. D 85, 055010 (2012) [arXiv: 1201. 2898 [hep-ph]]. 
[20] R. Nevzorov, arXiv:1205.5967 [hep-ph]. 

[21] D. J. Miller, A. P. Morais and P. N. Pandita, arXiv:1208.5906 [hep-ph]. 

[22] W. Porod, Comput. Phys. Commun. 153, 275 (2003) [arXiv:hep-ph/0301101]. 

[23] A. Djouadi, J. L. Kneur and G. Moultaka, Comput. Phys. Commun. 176, 426 
(2007) [arXiv:hep-ph/0211331]. 

[24] B. C. Allanach, Comput. Phys. Commun. 143, 305 (2002) [arXiv:hep- 
ph/0104145]. 

[25] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, arXiv:hep-ph/03 12045. 

[26] D. Chowdhury, R. Garani and S. KVempati, arXiv: 1109.3551 [hep-ph]. 

[27] G. Belanger, S. Kraml and A. Pukhov, Phys. Rev. D 72, 015003 (2005) 
[arXiv:hep-ph/0502079]. 

[28] G. Aad et al. [ATLAS Collaboration], arXiv:1207.7214 [hep-ex]. 

[29] S. Chatrchyan et al. [CMS Collaboration], [arXiv:1207.7235 [hep-ex]]. 

[30] S. F. King, S. Moretti and R. Nevzorov, Phys. Lett. B 650, 57 (2007) 
[arXiv:hep-ph/0701064]. 

[31] P. Athron, S. F. King, R. Luo, D. J. Miller, S. Moretti and R. Nevzorov, AIP 
Conf. Proc. 1200, 466 (2010) [arXiv:0909.4530 [hep-ph]]. 

[32] Y. Hosotani, Phys. Lett. B 126 (1983) 309. 

[33] J.L. Hewett, T.G. Rizzo, Phys. Rept. 183 (1989) 193; 

58 



[34] P. Langacker, Rev. Mod. Phys. 81, 1199 (2009) 

[35] J. Rich, D. Lloyd Owen and M. Spiro, Phys. Rept. 151, 239 (1987). 

[36] P. F. Smith, Contemp. Phys. 29, 159 (1988). 

[37] T. K. Hemmick et al, Phys. Rev. D 41, 2074 (1990). 

[38] B. A. Ovrut and H. J. Schnitzer, Nucl. Phys. B 184 (1981) 109. 

[39] L. J. Hall, Nucl. Phys. B 178 (1981) 75. 

[40] S. P. Martin and M. T. Vaughn, Phys. Lett. B 318 (1993) 331 [arXiv:hep- 
ph/9308222]. 

[41] K. Nakamura et al. [Particle Data Group], J. Phys. G 37 (2010) 075021. 

[42] H. Baer, J. Ferrandis, K. Melnikov and X. Tata, Phys. Rev. D 66 (2002) 074007 
[hep-ph/0207126]. 

[43] The ATLAS Collaboration, ATLAS-CONF-2012-033. 

[44] The ATLAS Collaboration, ATLAS-CONF-2012-037. 

[45] The ATLAS Collaboration, ATLAS-CONF-2012-041. 

[46] C. Rogan, arXiv:1006.2727 [hep-ph]. 

[47] The CMS Collaboration, CMS PAS SUS-12-005 (2012). 

[48] The CMS Collaboration, CMS PAS SUS-12-002 (2012). 

[49] The ATLAS Collaboration, ATLAS-CONF-2012-103. The ATLAS Collabo- 
ration, ATLAS-CONF-2012-104. The ATLAS Collaboration, ATLAS-CONF- 
2012-105. The ATLAS Collaboration, ATLAS-CONF-2012-109. 

[50] P. Athron, S. F. King, D. J. Miller, S. Moretti and R. Nevzorov, Phys. Rev. D 
84, 055006 (2011) [arXiv: 1102.4363 [hep-ph]]. 

[51] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 710 (2012) 49 
[arXiv:1202.1408 [hep-ex]]. 

[52] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710 (2012) 26 
arXiv: 1202. 1488 [hep-ex]. 



59 



[53] S. Chatrchyan et al. [CMS Collaboration], arXiv: 1206. 1849 [hep-ex]. 

[54] J. P. Hall and S. F. King, JHEP 0908, 088 (2009) [arXiv:0905.2696 [hep-ph]]. 

[55] J. P. Hall, S. F. King, R. Nevzorov, S. Pakvasa, M. Slier, R. Nevzorov, S. Pak- 
vasa and M. Sher, Phys. Rev. D 83, 075013 (2011) [arXiv:1012.5114 [hep-ph]]. 

[56] S. F. King and A. Merle, arXiv: 1205.0551 [hep-ph]. 

[57] J. P. Hall and S. F. King, JHEP 1106, 006 (2011) [arXiv: 1104.2259 [hep-ph]]. 

[58] M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. 
B 426, 269 (1994) [hep-ph/9402253]. 

[59] D. M. Pierce, J. A. Bagger, K. T. Matchev and R. -j. Zhang, Nucl. Phys. B 
491 (1997) 3 [hep-ph/9606211]. 



60 



